Dimensional transitions in small Yukawa clusters 



A. Radzvilavicius, 0. Rancova, and E. Anisimovas 
Department of Theoretical Physics, Vilnius University, 
Sauletekio al. 9-III, LT- 10222 Vilnius, Lithuania 
(Dated: March 3, 2013) 

We provide the detailed analysis ol structural transitions leading to the rapid changes in dimen- 
sionality ol small Yukawa clusters. These transformations are induced by the variations in the shape 
of confinement as well as the screening strength. We show, that even in the most primitive systems 
composed of only a few strongly interacting particles, the order parameter exhibits a power-law 
behavior in the vicinity of the critical point of the continuous transition. The critical exponent 
7 = 1/2 is found to be universal in all studied cases, which is consistent with the general theory of 
continuous phase transitions. 

PACS numbers: 61.46. Be, 64.60.-i, 52.27.Lw 



I. INTRODUCTION 

The emergence of ordered patterns within complex sys- 
tems of interacting entities draws attention of researchers 
in diverse fields, including physics, biology, mathematics 
and computer science [TH3]. Confined Wigner crystals 
are among the most primitive systems where the phe- 
nomenon of self-organization is observed. A small num- 
ber of charged particles is placed into a confining poten- 
tial well, where in the limit of low temperatures ordered 
structures are formed. This crystalline type of matter 
was successfully realized experimentally in a variety of 
well known systems, such as electrons on the surface of 
liquid helium |3J [5] , cooled ions in traps [B] or strongly 
coupled particle clusters in complex plasmas [Jj. 

Systems of strongly coupled particles are of high sci- 
entific interest due to various collective phenomena, e.g. 
cooperative dynamics, waves and phase transitions. First 
order transitions between solid and liquid states of mat- 
ter, namely melting and crystallization, were widely in- 
vestigated in the studies of many-particle dusty plasma 
crystals, also known as Yukawa clusters [HI [9] ■ Crystalline 
structures formed by dust particles in complex plasmas 
turned out to be an extremely handy tool for these stud- 
ies, as the convenient length and time scales, stability and 
transparency of these systems allow for direct optical ob- 
servation and accurate measurements [7J. On the other 
hand, there is another type of transitions, occurring in 
the simplest few-particle systems confined by asymmet- 
ric traps. These are observed when a small change in 
one of the control parameters causes a sudden change in 
the dimensionality of the system, and therefore are called 
dimensional or zigzag phase transitions [TU]- Note that 
these transitions take place in finite systems and are anal- 
ogous to the second order phase transitions commonly 
defined and studied in the thermodynamic limit. 

Dimensional transitions in small two-dimensional 
Yukawa clusters have been extensively studied in Ref. [TT] 
both experimentally and numerically. The authors 
demonstrated an excellent agreement between the com- 
puted configurations of particles and the arrangements 
observed in complex plasma experiments. Structural 



zigzag transitions were induced by variations in the num- 
ber of particles, the value of the Yukawa shielding pa- 
rameter k or the shape of the confinement potential. A 
power-law behavior of the order parameter in the vicin- 
ity of phase transition was observed, and interpreted as 
a characteristic feature of second order phase transitions. 
However, the numerical values of critical exponents pro- 
vided in pj] cast some doubt. These values are distinctly 
different from the classical mean-field value 1/2 which one 
would naturally expect in a finite system of just a few par- 
ticles. Therefore, one of the main motivations behind the 
current work is to provide the results of numerical mod- 
elling with higher precision and show the universality of 
critical exponents. We also extend the investigation to 
the clusters of different sizes and dimensional transitions 
of other types. 

Our paper is organized as follows. In section [TT] the 
model system is described and the procedure of our calcu- 
lations is presented. Section [llT| presents the results of the 
numerical modelling of dimensional transitions, grouped 
by their type in subsections. Main points of the article 
are summarized in section |IV| Additionally, Appendix 
[A| provides details on the analytically solvable cases and 
exact values of critical parameters. 

II. NUMERICAL PROCEDURE 

We investigate numerically systems of N identical par- 
ticles of mass m and charge Q, interacting through the 
Yukawa inter-particle potential. The interaction energy 
of two charges embedded in a screening environment thus 
reads 



where n stands for the shielding strength (inverse screen- 
ing length) and controls the range of interaction. Par- 
ticles are kept together by a harmonic confinement po- 
tential V c (r) = \muol (x 2 + a 2 y 2 + z 2 ) in 3D or V c (r) = 
\rnuJ 2 . (x 2 + a 2 y 2 ) in 2D. The anisotropy parameter a 
controls the shape of the confinement, which reduces to 
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the symmetric spherical or circular form at a = 1, and 
takes the shape of oblate (a > 1) or prolate (a < 1) 
spheroid in 3D and an ellipse in 2D. 

In the regime of strong correlations the potential en- 
ergy dominates over the kinetic one, and the total energy 
of the model system is given by 

N 1 N 1 

C7-(n, . . . , rjv) = J2 2 ( x i + a2 y* + Z D + Y, — e ~ Kr " ■ 

(2) 

The units of length and energy are conveniently chosen as 
r o = (Q 2 /47re ma;Q) 1 / 3 and E = Q 2 /47re r . Obviously, 
k is now measured in r^ 1 . In two dimensions, the z 
coordinates of all particles are set to zero, so that all the 
particles lie within the (xy) plane. 

Stable arrangements of particles correspond to the lo- 
cal and global minima of the potential energy ©. In 
our present work, stationary states are located by the 
method of multiple heating-relaxation cycles, based on 
the Monte Carlo and numerical minimization algorithms. 
The method was already proven to be efficient and reli- 
able in our previous studies [T5]. As it turned out, the po- 
tential energy landscape of ^ is rather complex even for 
small values of N and might be described as a collection 
of local minima, separated by potential barriers of various 
heights. In a first stage of the algorithm, thermalization 
takes place, that is, the system is heated to the temper- 
ature high enough to overcome all the potential barriers. 
This stage is accomplished by performing a large num- 
ber (few thousand) of Metropolis Monte Carlo steps [13] , 
which leads to a configuration, that may be regarded as 
drawn randomly from the Boltzmann distribution cor- 
responding to the given temperature. Each minimum 
controls a certain area of coordinate space, called basin 
of attraction. The area of attraction varies from basin 
to basin, which means that different stable states are re- 
alized with different probabilities [12j [14] . Some of the 
minima are located within the regions with steep walls, 
while others lie in broad shallow valleys and therefore re- 
quire a considerable effort to find. In a second stage of the 
computational procedure, the temperature of the system 
is suddenly set to zero and the closest local minimum of 
potential energy is located by employing the steepest de- 
scent and Newton optimization techniques. As frequently 
there is more than one local minimum [12] , the whole cy- 
cle of thermalization and relaxation is repeated a large 
number of times, to ensure that all of the basins are vis- 
ited and all stationary points of Q are revealed. 

Departure of the anisotropy parameter a from unity 
breaks spherical (or circular in 2D) symmetry and, as a 
result, prolate or flattened structures are formed. Even- 
tually, at the critical value of parameter a c , dimensional 
transitions are observed, three-dimensional clusters be- 
come planar, while two-dimensional structures are trans- 
formed into linear ones. In order to determine the critical 
values of a with high precision, we repeat our calcula- 
tions by incrementing a in small steps. As it will become 
evident shortly, properties of Yukawa clusters, including 



critical values of the anisotropy parameter and critical 
exponents, depend strongly on the screening parameter 
k. Therefore, in most cases we use four different values 
of Yukawa screening strength, k = 0, 1, 2, 3. In the sim- 
plest case of k — 0, inter-particle potential reduces to the 
simple unscreened Coulomb interaction. 

A second order phase transition is marked by a sud- 
den appearance or disappearance of some property of the 
system, called an order parameter, in response to a small 
change in a control parameter. We investigate dimen- 
sional phase transitions by keeping an eye on the total 
potential energy and order parameter (y), which is de- 
fined as the root mean square of the coordinate y: 

/ N vV» 

Naturally, dimensional transitions are signified by a sud- 
den change of (y) to zero. In particular, (y) is a good 
choice for an order parameter, since (y) = in the ID 
(2D) configuration and (y) > in the 2D (3D) config- 
uration. The state variables that determine the system 
configuration are then N, k and a. 



III. RESULTS 
A. 2D — > ID transitions 

We first investigate the simplest few-particle two- 
dimensional systems, undergoing 2D — > ID structural 
transitions. In the case of a = 1, the confinement poten- 
tial is symmetric and particles form ordered states, that 
were previously observed experimentally and modelled 
theoretically [T5J [TB]- Systems with N — 3, 4, 7 parti- 
cles form only one stable configuration (ground state), 
while clusters with N = 5, 6 particles in symmetric traps 
support both ground- and one metastable configurations. 
Various states can be represented by listing the occu- 
pation numbers of different shells — ground state of 5- 
particle system is therefore the configuration (0, 5) and 
metastable state is (1, 4) (for the arrangements of parti- 
cles see Figure[T]). As the anisotropy parameter a departs 
from unity, metastable states can become ground states, 
some states can disappear and new ones appear. In all 
investigated cases, however, there is only one stationary 
configuration near the dimensional transition — a zigzag 
shaped pattern, which soon becomes a ID linear chain of 
particles at a > a c . 

In the simplest symmetric case of N = 3, particles form 
an equilateral triangle in the (xy) plane. As the value of 
a increases, the triangular configuration is gradually de- 
formed until the transition occurs at a — a c ~ 1.55, as 
shown in Figure [I] with k = 0. In fact, dimensional tran- 
sition in three-particle system can be modelled analyti- 
cally, which gives the value of a c = ^jYljh (see Appendix 
|A[). Numerical simulation gives exactly the same value. 
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Figure 1. Value of the order parameter (y) as a function of 
the anisotropy a for 2D clusters with N = 3-5 particles. In- 
sets show typical arrangements of particles in various stages 
of compression. Note how quickly metastable state (1, 4) dis- 
appears (top left corner). 



The structural transition is slightly more intriguing in 
the case of N = 4 particles. As it can be deduced from 
the evolution of (y) in Figure [T] there is a discontinuity 
in a first derivative of the order parameter d (y) /da at 
the value of anisotropy parameter a ~ 1.69. As it turns 
out, there are two stages of the four-particle cluster com- 
pression. In the first, slow stage, four particles form a 
rhombus-shaped structure, with the particles located ex- 
actly on x or y axis. Later, the stage of rapid compres- 
sion takes over, with two particles departing from the line 
x = and forming a zigzag shaped pattern. Two typ- 
ical rhombus- and zigzag-shaped configurations are pre- 
sented in the insets of Figure [T] As it was shown in the 
previous study, the transition between these two stages 
is followed by the specific oscillation of the heat capac- 
ity [17] ■ Analogous scenario applies to the other clusters 
with even number of particles. The dimensional transi- 
tion is observed at a c « 2.04, where (y) suddenly drops 
to zero. 

Two competing configurations are first observed in 
the case of N — 5 particles, namely states (0, 5) and 
(1, 4). As it is depicted in Figure [I] metastable state 
(1, 4) exists only in the narrow window of anisotropy, 
1 < a < 1.05. On the other hand, the pentagonal ground 
state undergoes a continuous structural transformation, 
forms a zigzag-shaped cluster and finally becomes linear 
at a c « 2.50. 

Structural transitions become more complex for the 
systems with N > 6. Six particles in a symmetric con- 
finement can form two stable states. As a increases, 
metastable (0, 6) state vanishes at a = 1.05 only to reap- 
pear again and become a new ground state later. The 
former ground state (1, 5) then disappears completely 
near a = 1.22. Six- and eight-particle clusters both fea- 
ture the same discontinuity in d (y) / da as four-particle 



system, discussed above. 

We have already seen in Figure [l] that critical value of 
the parameter a c increases with N, when inter-particle 
interaction is of the Coulomb type. As it might be ex- 
pected, a c also grows as Yukawa potential screening pa- 
rameter k is increased, which is shown in Figure [2] for 
N = 3-6. The critical value of the anisotropy increases 
rapidly for k < 1.5 and almost saturates for high values 
of screening, i.e. k > 4.0. 

The lines in Figure [2] represent boundaries between 
different phases of clusters. The structures are two- 
dimensional below the line and form linear configurations 
above. Although we devote most of the present work to 
the transitions induced by deformations of the confine- 
ment well, structural changes can actually be caused by 
variations in any of three parameters a, k or N. As the 
figure shows, two-dimensional cluster can become linear 
without any changes in a, for example, when a value of 
k is diminished, or when a particle is removed from the 
system. 
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Figure 2. Phase diagram for two-dimensional Yukawa clusters 
of N particles. Below the corresponding line, the cluster of N 
particles is two-dimensional and one-dimensional above. We 
see, that dimensional transitions can be induced by variations 
in any of three parameters a, k or N. 

We further examine the power-law behavior of the or- 
der parameter (y) near its critical point a c in more detail. 
The power law is easily identified by plotting the loga- 
rithm of the order parameter, lg((y)), as a function of 
lg(a c — a). The function turns out to be linear for small 
values of (a c — a). This observation confirms that in the 
vicinity of the transition point, the order parameter (y) 
demonstrates a power-law behavior, which is a typical 
property of second-order phase transitions: 

(y) cx {a c - a y. (4) 

We determine the values of the exponent 7 near the 
critical point by analyzing the slope of the above dis- 
cussed log- log plot. Namely, we take the numerical 
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derivative of the function lg ((y)) = / (lg (a c — a)). Cal- 
culated exactly at the critical point this derivative yields 
the exact 'theoretical' value of the critical exponent. 
However, in an experimental or numerical investigation 
the precise location of the critical point may not be 
known. Thus, calculating the numerical derivative a bit 
away from the critical point we are able to mimic the 
uncertainty and errors present in a realistic experimental 
situation. It turns out, that in all cases, 7 = 1/2 as long 
as a is close to its critical value a c (Figure [3]). However, 
the local value of the exponent (determined as the nu- 
merical derivative) is very sensitive to the deviation of 
anisotropy parameter from a c . 

Figure [3] shows the dependence of the power-law expo- 
nent 7 on the deviation of a from its critical value. The 
case of 2D cluster with N — 3 particles and four differ- 
ent values of screening length is presented. We see, that 
7 departs from the value of 1/2 significantly when the 
deviation from a c reaches third decimal and attains its 
minimum near the first decimal. Furthermore, the expo- 
nent 7 attains significantly lower values far from a c in the 
systems with stronger screening. Other than that, there 
are no qualitative differences in the critical behavior of 
systems with different values of k. The exponent 7 of the 
other systems with N > 3 behaves similarly to the case 
of three particles presented here. 
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Figure 3. Critical exponent of the transition observed in 2D 
three-particle systems as a function of deviation of a from its 
critical value a c for different strengths of screening. 

A thorough analysis of transitions in a planar cluster of 
five particles was reported by Sheridan in . The given 
value of critical parameter is a c = 2.96 and the critical 
exponent is said to be 7 = 0.39 7^ 1/2, while we find 
critical anisotropy parameter to be a c = 3.01. According 
to the results of our calculations, a value of a = 2.96 
corresponds to the exponent 7 = 0.37, which is close 
to the value reported by Sheridan et. al. Therefore, the 
reason of differences between the results, published in [TT] 
and our discoveries almost undeniably lies in the extreme 



sensitivity of (y) on the value of anisotropy parameter 
and distance from its critical value. Thus, the accuracy 
of the results presented in [TT] is probably not sufficient. 

As it was demonstrated in [llj and shown in Figure 
[2] continuous transitions might also be induced by the 
variations in the screening strength k. By keeping N 
and a constant, we gradually change a value of k while 
tracking the changes parameter (y) undergoes. The di- 
mensional transition takes place when the value of (y) 
suddenly drops to zero, at which point the critical value 
of k is obtained (see the inset of Figure [4]). It turns out, 
that (y) exhibits the same power law behavior near the 
transition, i.e. (y) oc (n — k c )P. Figure [i] shows the de- 
pendence of critical exponent /3 on a logarithm of the 
distance from the critical value k c for three systems. As 
opposed to the results presented in [11], we see again, 
that in all cases close to the transition point f3 = 1/2. 
Moving away from the critical point, however, exponent 
(3 departs from the value of 1/2 significantly. 
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Figure 4. Critical exponent /3, as a function of deviation of k 
from its critical value k c for three 2D clusters. The evolution 
of the order parameter (y) close to the dimensional transition 
(inset). 



B. 3D -> 2D transitions 

We further investigate structural transitions in three- 
dimensional Yukawa clusters with N — 4 to N — 8 par- 
ticles and integer values of screening parameter up to 
k = 3. Increased values of the parameter a, turn initially 
spherical structure into the oblate one; eventually, after 
the anisotropy parameter reaches its critical value a c , di- 
mensional phase transition takes place and familiar two- 
dimensional clusters are formed. In three-dimensional 
transformations of five- and six-particle clusters, two dif- 
ferent final states are possible, as opposed to the zigzag 
transitions in 2D, where only one linear configuration can 
be formed. Therefore, in 3D — > 2D transitions, there is 
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a distinct value of a c for each final configuration and we 
are concerned by the properties of phase transitions of a 
particular stable state. 

As the confinement potential well is squeezed in y di- 
rection, it is handy to label small clusters according to the 
arrangement of particles in the projection to (xz) plane, 
in a manner similar to the state labeling by shell occupa- 
tion numbers in two dimensions. Moreover, particles in 
three-dimensional anisotropic traps frequently organize 
themselves within the layers, parallel to the (xz) plane. 
For the sake of clarity and unambiguous definition of the 
configurations, we will also use a list of particle numbers 
in distinct layers, enclosed within curly brackets. 

The simplest system, undergoing a non-trivial 3D — > 
2D transition is the cluster composed of four particles. 
Not surprisingly, four particles in a symmetric three- 
dimensional trap form a regular tetrahedron, and there 
is only one possible square-shaped (0, 4) state in two di- 
mensions. Figure [5] shows dependence of the order pa- 
rameter (y) and potential energy of a system E on the 
anisotropy parameter a. We see, that (y) changes con- 
tinuously and the transition is remarkably similar to the 
one in 2D case of N = 3 particles. The potential energy 
gradually increases as the potential trap is flattened, un- 
til a two-dimensional structure is formed at a c « 1.22. 
As it is demonstrated in Appendix [XJ this symmetric 
transition can be modelled analytically; the critical value 

1 /2 

turns out to be a c = (4V2/(1 + 2V2)) — exactly the 
same as determined in our numerical modelling. Natu- 
rally, the value a c is sensitive to the range of the inter- 
particle Yukawa potential. As Figure [H] shows, critical 
value of the anisotropy parameter increases rapidly with 
the strength of screening for k < 2 and significantly 
slower after that, thus reminding of the transitions from 
two- to one-dimensional configurations. 
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Figure 5. Order parameter (y) and potential energy £ of a 
four-particle 3D Coulomb cluster in the asymmetric potential 
trap with anisotropy parameter a. 

As it was already pointed out before, there are two 



competing stable states observed in a two-dimensional 
system with N = 5 particles. Five particles in a spheri- 
cally symmetric 3D confinement potential, however, can 
form only one stable configuration. Slightly increased 
anisotropy leads to the formation of a three-layer struc- 
ture, with the arrangement of particles within these lay- 
ers being {2, 2, 1}. As parameter a increases above the 
value of a = 1.05, two layers merge forming a square 
and thus transforming the configuration into a pyrami- 
dal structure {4, 1} with projection (1, 4) xz . This struc- 
tural transition from three-layered cluster to the pyrami- 
dal configuration is signified by the discontinuity in the 
derivative d (y) /da (Figure [6| . 

As it is demonstrated in Figure [6] for a pure Coulomb 
interaction, a second stable state appears when the 
anisotropy parameter reaches the value of «o ~ 1.29. 
A new pentagonal state undergoes an asymmetric di- 
mensional transition and is soon transformed into the 
new ground state (0, 5). The metastable configuration 
(1, 4) xz , on the other hand, becomes two-dimensional 
only at a c s» 1.60, through the so-called "pyramidal" 
transition mechanism. 

Both point of appearance of the second state in five- 
particle system ao, and its critical value a c depends on 
the type of the interaction potential and its screening 
parameter n. As it is shown in Figure [HJ both parame- 
ters grow with the strength of screening. The distance 
between ao and a c , however, rapidly diminishes. As 
the screening reaches the value of K ~ 4.5, two lines 
merge and a new stable state appears already in its two- 
dimensional pentagonal form (0, 5). 
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Figure 6. Order parameter (y) and potential energy E of a 5- 
particle 3D Coulomb cluster in the asymmetric potential trap 
with anisotropy parameter a. Here and later the insets show 
projections of configurations to (xz) plane. 

A pyramidal configuration might be described as a pla- 
nar base, composed of n = 4-6 particles lying parallel to 
the (xz) plane, and a single particle located right above 
the center of the base, that is, configuration {N — 1, 1}, 
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(1, N — \) xz . A pyramidal structural transition takes 
place, when the apex of the polyhedron is pushed into 
the base, thus becoming a two-dimensional configuration 
with only one particle in the center. A typical behavior 
of the order parameter (y) during such transitions was 
already discussed and is presented in figures [6] and [7j As 
a matter of fact, the dimensional transitions of a pyrami- 
dal type can be modelled analytically and exact values of 
critical parameters a c can be found, see Appendix |A"| 

Even more stable configurations are observed in clus- 
ters with N = 6 particles, as Figure [7] shows for the 
Coulomb inter-particle potential. Evolution of the sys- 
tem starts with a single stable state in the symmetric 3D 
trap — the octahedral configuration (full line in Figure 
[7b. As the parameter a increases, this bipyramid is de- 
formed by pushing two of its particles lying exactly on 
y-axis, towards each other, thus lowering the height and 
forming a configuration {1, 4, 1} with projection (1, A) xz . 
(y) decreases slowly, until the said two particles start to 
depart form the y-axis near a rs 1.46, at which point a 
phase of rapid deformation begins. Unfortunately, right 
after this happens, the stable state disappears. The same 
scenario of bipyramidal deformation also applies to the 
larger clusters, e.g. N — 7, 8, 9, and has a specific, 
well recognizable shape of its (y) = f(a) curve, with the 
segments of slow and rapid changes. 
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Figure 7. Order parameter (y) and potential energy E of the 
six-particle 3D Coulomb cluster in the asymmetric potential 
trap with anisotropy parameter a. 

A new metastable state emerges near a w 1.06: the 
particles lie on the six vertices of two parallel equilat- 
eral triangles, centered precisely on y-axis and rotated 
by 7r/3 with respect to each other, that is, state {3, 3} 
and (0, 6) X z- These triangular layers are pushed towards 
each other by deformations of the confinement, however, 
they fail to ever become a truly two-dimensional con- 
figuration. Instead, as it is demonstrated in Figure [7J 
the configuration ceases to exist at a 1.52 where the 
r.m.s. value of y coordinate is still (y) ~ 0.05 > 0. How- 
ever, right before the disappearance, a new similar purely 



two-dimensional state shows up. The new planar config- 
uration is composed of six particles lying on the vertices 
of two triangles of slightly different sizes (see Figure [7]) . 
Therefore, in a brief range of a values these two states 
exist simultaneously and there is no continuous transi- 
tion between them. Finally, a pyramidal configuration 
{5, 1} appears near a ss 1.19 and undergoes the usual 
pyramidal dimensional transition at a c ~ 1.59, the value 
predicted by our analytical model (appendix [A|. 
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Figure 8. Critical value of anisotropy parameter a c as a func- 
tion of re in three-dimensional configuration with N = 4 par- 
ticles and two states of N = 5 clusters. The line marked as 
«o corresponds to the appearance of the metastable state in 
a five-particle cluster. 

Close to the critical point of continuous transitions 
from three- to two-dimensional systems, a power-law be- 
havior of the order parameter is detected once again, i.e. 
(y) oc (a c — a) 1 . In the same manner as in 2D case, 
Figure [9] shows the dependence of the power-law expo- 
nent 7 on the logarithm of a c — a. It turns out, that in 
a close vicinity of transition point, the critical exponent 
7 = 1/2 does not depend on the screening strength n. 
Deviations from this value occur when the departure of 
a from its critical value reaches third decimal. Just as 
in the two-dimensional case, the value of a critical ex- 
ponent is lower for systems with stronger inter-particle 
potential screening, and drops as low as 7 » 0.35 for 
k = 3. Essentially the same behavior of the exponent 
7 is observed in larger three-dimensional systems, where 
dimensional transitions take place, be it pyramidal tran- 
sitions or transformations of any other type. 

As the number of particles grows, more and more stable 
states emerge and, as a consequence, (y) — f(a) graphs 
become convoluted and somewhat difficult to study. An 
illustrative example is given in the inset of Figure |10[ 
where the behavior of the order parameter in stable states 
of 20-particle Coulomb system is presented. We can 
still see a few distinct continuous phase transitions in 
the vicinity of a w 2.26, however different lines become 
hardly distinguishable at the lower values of a. In very 
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on the root mean square of the distance from y-axis: 




Figure 9. Critical exponent 7 for a dimensional transition, 
observed in 3D four-particle systems as a function of devia- 
tion of a from its critical value a c for different strengths of 
screening. 



large systems, the values of (y) for all metastable states 
lie virtually on the same line, as it is illustrated in Figure 
10 with a 100-particlc cluster and K = 0. 
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Figure 10. Order parameter (y) of large Coulomb systems 
with N = 100 and N = 20 (inset) particles as a function of 
the anisotropy a. 
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As Figure [TT] shows for N = 3-6, dependencies of 
the order parameter (p) on the anisotropy a are not 
smooth and in some cases feature discontinuities. We 
conclude, that in fact there is no direct transition from 
three- to one-dimensional configurations. Instead, the 
system is first transformed into the elongated 2D zigzag 
pattern, and only later, 2D — > ID structural transition 
takes place. With that in mind, there is no surprise, that 
the values of critical parameters a c found by lowering a 
are the exact inverses of those, determined in subsection 

MM 




Figure 1 1 . Order parameter (p) as a function of the anisotropy 
a < 1 for Coulomb clusters with N = 3-6 particles. 

Large three-dimensional clusters in prolate traps be- 
come one-dimensional through the mechanism, which 
seems to be universal for all values of N used in our 
modelling. At first, the system is squeezed and elongated 
until the particles arrange themselves into the shape of 
double-helix. As a is lowered further, a number of heli- 
cal turns decreases, until the helix unwinds and cluster 
becomes a two-dimensional zigzag configuration. 2D sys- 
tem then undergoes the usual zigzag transition with the 
power-law behavior near the critical point. 



It might be worth discussing the structural evolution of 
Yukawa clusters confined by traps with a prolate equipo- 
tcntial surface. In our model, this effect is achieved by 
lowering the value of anisotropy parameter a towards 
zero. In that way, elongated clusters are formed, with low 
potential energies and high values of (y). Consequently, 
in order to study this type of structural transitions, a 
new order parameter must be defined. We choose to rely 



IV. CONCLUSION 

Confined Yukawa clusters are among the physical sys- 
tems, where simple interparticle interactions lead to 
the emergence of complicated patterns and spontaneous 
ordering. In this article, we present our findings in 
the numerical and analytical studies of two- and three- 
dimensional clusters confined by asymmetric parabolic 
traps. 
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We confirm, that dimensional transitions from oblate 
three- to two-dimensional systems as well as from planar 
to linear configurations can be induced by changes in the 
anisotropy of the confinement a, and screening strength 
k. On the other hand, there are no direct transitions from 
three- to one-dimensional systems in prolate harmonic 
traps; two-stage transformations take place instead. 

A critical value of the anisotropy parameter in general 
grows with the screening strength k. The growth is steep- 
est for small values of k and almost saturates for the large 
ones. In a close vicinity of dimensional phase transition, 
the order parameter (y) exhibits power-law dependence 
on a control parameter, be it a or n. 

In all cases studied here, the critical exponent is found 
to be universal and equal to 1/2, which is consistent 
with the general theory of second order phase transitions. 
However, a value of the power-law exponent turns out to 
be very sensitive to the deviations of a control parameter 
from its critical value. Far from the critical point, the 
exponent attains lower values in systems with stronger 
screening and shorter range of the inter-particle interac- 
tion. 



Appendix A: Analytical values of a c 

In a few simplest cases of high symmetry, when total 
energy of a cluster after a transition depends on a single 
generalized coordinate, values of a critical anisotropy pa- 
rameter a c can be identified analytically. These are the 
transitions in 3-particle 2D system, 4-particle 3D cluster 
and all of the pyramidal transitions. 

Consider the iV-particle three-dimensional system un- 
dergoing a pyramidal 3D — > 2D dimensional transition. 
Immediately after the transition, a planar cluster consists 
of n = N — 1 particles, positioned on the circumference of 
a circle with radius R, and a single particle in the center 
of confinement. The 2D cluster lies in (xz) plane. The 
total potential energy of the system can be expressed as 



U (R) 



1 



I? 



/(") 



n 

R' 



(Al) 



A second term here represents the Coulomb interaction 
energy of n particles positioned on the circle with radius 
R, thus forming a regular polygon. The function f(n) 
depends only on the number of particles and is 



1 1 

4 + 2 



a/2-1 

sin -1 (rror/n) if n is even, 



/(») = 



(n-l)/2 



- sin 1 (rrnr/n) if n is odd. 



m— 1 



,(A2) 

By setting 8Uq jdR = we find the equilibrium radius to 
be simply 



i?=(l + /(n)) 1/3 . 



We further perturb the system by shifting the base by 
the small distance of Sy, and, keeping the center of mass 
at the center of the confinement, an apex by the distance 
of nSy in an opposite direction. Now, the total energy of 
the system is 

U v = \n{R 2 + a 2 5y 2 )+ l -a 2 {n8y) 2 

n , . n (A4) 

H f(n) H . 

R y/(R 2 + (n+l) 2 5y 2 

By expanding the right hand side in powers of Sy up to 
the quadratic term, we get 



17, 



1 



P ~ U a + -not Sy (1 + n) 



n(n + l) 2 6y 2 
2R? 



By requiring that 

xtt 1 2r 2/ n , \ n(n + l) 2 5y 2 
SU = -najy (l + n) — 







(A5) 



(A6) 



for any small Sy, and making use of equation |A3| we find 
the critical value of anisotropy parameter to be 



N 



f(N-l) + l 



(A7) 



Table I. Values of the function / and critical parameter a c 
of 3D Coulomb systems undergoing dimensional transition of 
the pyramidal type. 



N 


f(N-l) 




5 


i + f 


2 J w 1.5983715 

V 5 + 2V2 


6 


H 


l/3\/5 + 2\/5 - V5 « 1.5889766 


7 


5 1 

4 + 7I 


2 J ^— w 1.5734727 

Y 27 + 4^ 


8 


2.3047649 


1.5558750 



(A3) 



Values of the function f(n) and corresponding critical 
parameters of pyramidal transitions are collected in table 
[I] for all 3D clusters with this type of structural transfor- 
mation. We see, that in general, a c slightly decreases 
with N in the range of 5-8 particles. 

In the two-dimensional case of N = 3 particles, the 
system forms a triangular cluster. During the structural 
transition, one of its particles is pushed in-between the 
others, thus forming a linear structure. This case is ba- 
sically a generalization of pyramidal transitions to the 
two-dimensions. Therefore, equation |A7| is still valid and 
we find critical anisotropy to be a c = W12/5, which is 
exactly the value observed in our numerical modelling. 

A slightly different mechanism of transformation is ob- 
served in 3D — > 2D transition of highly symmetrical 
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N = A particle cluster. The final configuration is a 2D 
square, with potential energy 



U (R) 



1 



-NR Z 



N 
R 



W). 



(A8) 



Minimization of A8 by solving dUo /OR = in turn gives 
R = f(N) 1 ' 3 . Right before the transition, two particles 
sharing a common diagonal in the final square are ele- 
vated by the distance of Sy above {xz) plane, while other 
two are located at the same distance below it; ergo the 
energy of the perturbed system 



Ur, 



N 



y/2R 2 + 48y 2 



Un 



2 y 2i? 3 



(A9) 



By setting SU = U p — Uq = 0, we again get 




1.2155625. 



(A10) 



This value is exactly the same as found by our numerical 
procedure. 
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